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' If massive neutrinos possess magnetic moments, they can undergo 

i | spin flip in a magnetic field. The magnetic fields needed for a 

meaningful measurement of neutrino moments could be very high 
and may occur in astronomical objects such as some supernovae 
or active galactic nuclei: they are typically chaotic ones. We de- 
velop the general theory of the passage of neutrinos through such 
fields. We also develop a simple model which becomes solvable 
in the high energy limit. Both helicities occur with equal prob- 
ability, independently of the initial distribution. Observational 
consequences are discussed. 
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1 Introduction 

The existence and magnitude of magnetic moments of neutrinos is one of 
the important outstanding questions of particle physics. It is well known 
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that in the standard model of electroweak interactions neutrinos are purely 
left handed; hence, they are massless and cannot possess magnetic moments 
either. However, in some minimal extensions of the standard model and also 
in in some grand unified models, e.g. in the SO(10) model, there are both left 
and right handed neutrinos present. Therefore, there is the possibility giving 
them mass and also, chirality flipping (Pauli type) interactions with the 
electromagnetic field. In such models, higher order electroweak interactions 
give rise to anomalous magnetic moments via loops in which a W and a 
charged lepton is circulating. This was first noticed by Lee and Shrock |l[ 
a long time ago. The magnetic moment thus arising, however, is a very 
small one, due, mainly, to the large mass of the W boson. It is given by the 
expression, 

3G F em u 19 m v 

a u = ■= ps 3 x 10 — — — . 

8V2tt 2 leV/i B 

Here Gp stands for the Fermi coupling constant and hb is the Bohr magneton. 

There exist fairly reliable upper limits on the masses of neutrinos; using 
those, the magnitudes of the moments predicted by this formula come out to 
be several orders of magnitude lower than the current observational upper 
limits, see, e.g. |2| for a recent review. (Current limits are typically around 
10~ 12 /is for flavor diagonal moments and a few orders of magnitude higher 
for transition moments.) As a consequence, any observation of a neutrino 
magnetic moment indicates the presence of some significant physics beyond 
the standard model. Using the observational upper limits on the moments, 
however,one concludes that such a measurement is likely to take place by 
means of a neutrino telescope. Indeed, the characteristic length associated 
with a magnetic dipole moment of magnitude fi placed in a magnetic field of 
strength B is given by L = 1/fiB. (Roughly speaking, this is the distance 
over which helicity flip occurs with a substantial probability.) In convenient 
units, 

- — ^ 3 x 10 . 

km fiB 

For a magnetic moment of 10 -12 /zs moving even in a field of 1 Megagauss, the 
characteristic length is about 3 x 10 4 km. Clearly, no terrestrial experiment 
can be designed for such a measurement. By contrast, an active galactic 
nucleus (AGN) or a supernova, to quote but two examples, provides us with 
the right laboratory for measuring small magnetic moments either because 
of the high magnetic fields and/or by the large size of the object in question. 

The trouble is that magnetic fields in astronomical objects are generally 
chaotic ones and the description of the motion of a magnetized neutrino in 
such an environment is hard to describe by the methods particle physicists 
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are used to. The main purpose of the present work is to develop the general 
formalism necessary for solving such a problem. Section [| is devoted to 
the development of this formalism, while in the subsequent one (Sec. |3|) we 
describe the relevant dynamics in front form. A solvable model (involving a 
reduced number of neutrino flavors) is described in Section f|. Finally, the 
results and their observational consequences are discussed in Sec. || 

2 General Formalism 

The problem of investigating the behavior of a microscopic system in a ran- 
dom environment is often faced in condensed matter physics. The classic 
example is, of course, Brownian motion, but similar problems occur in the 
statistical theory of turbulence and many other situations as well. Typically, 
such systems are modeled by an equation of motion of the form: 

d t X + F[X] = f, (1) 

where X is an element of a vector space (the phase space of dynamical 
variables), the functional F[.] is often, but not necessarily, local and time 
independent. Finally, / represents a random force perturbing the otherwise 
deterministic system. 

The modern theory of such systems was developed by Martin, Siggia and 
Rose || using a canonical formalism and by De Dominicis and Peliti Q 
using a functional integral formalism. For the purpose of generalization to a 
relativistic system, the formalism developed in ref. [|J is more suitable. 

There is an important difference, however, between systems described by 
an equation like ([I]) and ones we are interested in. In ([!]) the random 
perturbing force appears as an inhomogeneous term on the right hand side 
of the otherwise deterministic dynamical equation, whereas for a particle in 
a random magnetic field, the latter appears typically in the form B • Syj, 
where S is a spin operator and (p is an amplitude describing the propagation 
of the particle^. Hence, the formalism has to be modified in order to take 
this circumstance into account. 

We use a Hamiltonian formalism; manifest covariance is of little value, 
since the rest frame of the astronomical object provides us with a preferred 
frame of reference. The Hamiltonian is taken to be of the form: 

H = H + fiM ■ B, (2) 

2 In this work we are interested in neutral particles, hence interaction terms proportional 
to the electric charge are absent. 
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where Hq governs the propagation of the particle in the absence of a magnetic 
field. The quantity M describes the coupling of the particle to the magnetic 
field and it may contain degrees of freedom other than just spin: e.g. various 
flavors a particle comes in may have different types of magnetic couplings. 
Correspondingly, ip is a vector in the space of helicities and flavors. Finally, 
fi is the magnitude of some average of the magnetic moments involved: it 
carries the necessary dimension and the absence of a magnetic coupling can 
be described as the limit fi — > 0. 

The equation governing the propagation of the particle thus reads: 

d t (p + H <p + fiM ■ By? = (3) 

This equation is a deterministic one: there is no randomness involved. As- 
suming that one finds a solution of (0) for an arbitrary magnetic field and 
one is given the distribution of B over the ensemble of the magnetic fields 
of interest, one can determine the average behavior of the solution as well as 
the fluctuations around the average. 

Formally, let if (x|B) denote the solution of (H) as a functional of the 
magnetic field, where x stands for the spatial coordinates and the time. 
Assume further that one is given the characteristic functional, Z [j] of the 
distribution of the magnetic fields over the ensemble. Then the average 
behavior of (p is given by the expression: 

(<p{x\B))=<p[x\^jz\i] j=0 . (4) 

The fluctuations around the average are given by similar expressions involv- 
ing higher functional derivatives. 

Equation (|) can be useful provided two things are specified: 

• The characteristic functional, Z [j] 

• The solution of equation (|3]) for an "arbitrary" magnetic field. 

It is clear that, for any system but the very simplest ones, the two items 
specified above cannot be given in a closed form. (For a simple system to be 
discussed in Section |], one can give analytic expressions in the limit of high 
energies.) What one needs instead is a meaningful approximation scheme 
permitting the calculation of the quantities needed. 

First of all, we notice that for any physically acceptable ensemble of 
magnetic fields, the characteristic functional can be written as a functional 
integral, viz. 

Z[j] = j VB exp —S [B] x exp z(j • B). (5) 
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Here S stands for a generalized entropy functional and (...) denotes inte- 
gration over continuous variables and summation over discrete ones. The 
functional measure is normalized in such a way that Z [0] = 1. The general- 
ized entropy has two important properties: 

1. S> 

2. ft£Z[j] = 0. 

The first property follows from the fact that Z is the characteristic functional 
of a probability distribution, the second one means that the magnetic field is 
solenoidal, as it should be. In general, the entropy functional is determined 
from solving Maxwell's equations for some random current distribution; in 
the next section we shall consider a simple case of a static field. There- 
after, one has to decide on physical grounds the expression of the entropy 
functional. 

The simplest, physically reasonable approximation is to consider a Gaus- 
sian distribution, 

S=±(BC- 1 B), (6) 

where C is the correlation operator. In what follows, we shall assume a 
distribution of this type. 

Next, we construct the generating functional for the moments of p. We 
have to average over the manifold of solutions of equation so we have: 

Z[u\ = J VBV(pexp-^(BC- 1 B)expi(u<p) 

x 5 [id t p + H cp + /iB • Mcp] Det (id t + H + //B • M) , (7) 

where u is the source of p. The determinant is inserted in order to get 
the correct measure on the manifold of solutions. One can use the integral 
representation of the ^-functional, 5[.] in order to cast (0) into a form to 
which techniques familiar from quantum field theory are directly applicable. 
We write: 

5 [id t (p + H p + fiB ■ Mtp] = J Vip ] exp z(y? f (id t cp + H Lp + /iB ■ My?)) (8) 

It is to be noted that in tp^ is not the Hermitean conjugate of ip, but 
an independent functional argument; the notation is used because of the 
transformation properties of p^ under any symmetry group the argument of 
the delta- functional may have. Using (||) the generating functional for the 
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moments of ip can be brought to a form familiar from quantum field theory, 
viz. 



u, «' 



= JvBVifVif^exp-^BC^B) 

x expi((if^u + wV)) 

x expi{if^ (id t ip + H ip + /xB ■ M_ip)) 

x Det (id t + H + fiB ■ M) . (9) 

(We introduced a source for if* as well, in order to be able to generate Green's 
functions.) 

Equation (§) is recognized as being formally equivalent to a quantum 
field theory for the variables ip and ip\ the variable B merely serving as 
a mediator of interactions. The functional determinant can be handled by 
familiar techniques, typically, by introducing Fadeev-Popov ghost fields. The 
Gaussian integration over B can be carried out explicitly; however, this may 
not be the most desirable form of @. 

We conclude that the problem of describing the propagation of a neutral 
particle (a neutrino in particular) in a random magnetic field can be handled 
in a way familiar in quantum field theory. Familiar approximation methods, 
such as loop expansions, etc. are readily applicable to this problem as they 
are to any quantum field theory. 

So far, we have not specified the form of the dynamics entering the pre- 
ceding equations. This is an important question: casting the dynamics into 
an appropriate form leads to significant simplifications. 



3 Dynamics in Front Form. 

We use the front form of dynamics ||. As explained in a previous work 0, 
this formulation of dynamics is advantageous in a situation in which one 
considers the propagation of high energy particles (E 3> m, where m is the 
rest mass) and in which certain discrete symmetries, such as C and P play 
no significant role. Clearly, the propagation of high energy neutrinos falls 
into this category. In what follows, we assume an arbitrary flavor structure: 
all physical observables (masses, magnetic moments, etc. ) are matrices in 
flavor space. However, flavor indices need not be explicitly exhibited. 

We begin with the usual Dirac Lagrangian of a particle in an external 
electromagnetic field, F^ v : 

L = i> (V^ + m + ^F""<v) ^ (10) 
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We work in the rest frame of the magnetic field. Assuming the field to be 
a static oneQ, we can set F 0i = 0, Fy = €ij k B k . In the case of interest one 
has to solve the Dirac equation in an arbitrary static magnetic field, since we 
want to average the solution over an ensemble of the Bi. No explicit solution 
is known for such a problem. However, we proceed to show that in the high 
energy limit the problem can be solved in a closed form. 

We introduce a coordinate system in which two of the coordinates are null 
directions corresponding to characteristic lines of a relativistic wave equation, 
viz. : 



t = —= (x° -x 3 ) ,z = —= (x° + x 3 ) and x A ; (A = 1,2). 



'Ill 



v/2 V r V2 

Correspondingly, the metric is of the form, 

9zt = 9tz = 1, 9ab = Sab, (12) 

and all other components vanish. The momentum components conjugate to 
z and t are, 

-\- ~^ 

k~p 3 V2, ^ 3 ~* 00 ' ^ finite )' ( 13 ) 

respectively. In this equation, p is the momentum transverse to p^. 



A Dirac spinor can be decomposed along the null directions given in (11 
by introducing the mutually orthogonal projectors, 

Pt = \ltl\ P z = \lzl z (14) 
In what follows, we use the shorthand, 

= Prf, X = Pzi> (15) 

It is a straightforward matter to decompose ( |T0D according to the con- 
jugate null directions and express it in terms of the variables and x- The 
purpose of such an exercise is a very simple one. If, for the sake of definite- 
ness, t is regarded the "time" variable describing the dynamics of the system, 
only cf) obeys an equation containing d t . Hence, the conjugate component of 
the Dirac spinor obeys only an equation of constraint. The constraint can 
be, in turn, solved before one attempts to attack the problem of dynamics. 

3 In physical terms, this means that the characteristic time scale of change of the field 
is large compared to the time of passage of neutrinos. It is a straight forward matter to 
generalize the formalism for arbitrary electromagnetic fields. 
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After carrying out the decomposition of fllCf ) according to the null direc- 
tions, one finds: 

L = y/2[tf (id t -iV2fMe AB j A B B 

+ X ] fa - iV2fie AB lA B B ) X 
1 

+ 73 



lV (ij A d A + m - 4|/^3 £ab1 A 1 B j X 



(16) 



Variation of (|1^) with respect to x gives the constraint: 

(iy/2d s - 2tfie AB lA B B ) X + ^7* (V^ + ^e ABl A 1 B B 3 + m^j = 0. 

(17) 

The solution of flT7D with the correct Hermiticity properties, c/. 0] reads: 

- | cfe'<z^ |z ') -^ 7 * (i<y A d A +m- ¥^=B 3 e^ 7 V) CO 

(18) 

Here all arguments in which fll8|) is local, (i.e. x" 4 and t) have been sup- 
pressed. The matrix element of the operator Q is given by: 

(z\Q(B A ) \z>) = -j=exp(^V2 £dz"e ABl A B B (z")^ \e{z-z') (19) 

Let us notice that solving the constraint eliminates two components of 
the original, four component Dirac spinor. Therefore, instead of the original 
Dirac matrices one can go over to 2 x 2 Pauli matrices. One easily verifies 
that —ie AB ^ B — > o A gives the correct representation. We also introduce the 
Hermitean operator, p A = —id A for the transverse degrees of freedom. 

It is now a matter of straightforward algebra to use fll8|) in order to 
eliminate the variable x from the Lagrangian. We merely quote the result; 
it can be conveniently written in Hamiltonian form as follows. 

L = ird t (j)-H 
H = -2/j,^a A B A (j) 

+ f (-ia B e BC p c + m - /iV2B 3 a 3 

x Q(-B A 



(20) 



x (ia R e RS p s + m - fj,V2B 3 a 3 ^ 
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The canonical momentum is given byvr = iy/2tf. (Of course, the odd looking 
factor of y/2 in the definition of the canonical momentum can be eliminated 
by rescaling the time variable.) 

In ( p0|) we omitted all the symbols of integration over the longitudinal 
coordinate, z; terms in which the operator Q does not appear explicitly are 
local in all variables. 



The Hamiltonian appearing in fl20l) is exact. However, it is given by a 



rather complicated, non local and non linear expression: this is the cost we 
have to pay for explicitly eliminating the constraint. We now argue that one 
can introduce physically reasonable simplifications, as a result of which the 
problem becomes a manageable one. We notice that the exponential appear- 
ing in (^), regarded as a matrix in the space of spinors is unitary, hence all 
matrix elements are bounded. Therefore, one expects that at large values of 
| z — z'\ the exponent oscillates rapidly and contributes little to the Hamil- 
tonian. The dominant contribution is, probably, coming from small values 
of the difference of longitudinal coordinates. Hence, it appears reasonable to 



approximate the exponential in (19) by 1. It is to be emphasized, however 
that the linearization of the effective Hamiltonian in the magnetic field is 
not necessarily a harmless approximation. Many instances are known when 
a similar truncation of a dynamical system leads to the loss of important 
physics. The validity of this approximation therefore needs further study. 

We now observe that, after having made the approximation just discussed, 
there remain three types of terms in (ffi|): 

1. The local term, cx o a Ba- this is purely a helicity flip term in terms of 
the lightlike helicity defined by o~ 3 . 

2. Transverse spin-orbit coupling terms , oc <ta£ AB Pb'- these terms also 
give rise to the transverse kinetic energy term, since ct A €abP B vr£ RS Ps — 
PaP A - 

3. "Mass" terms, oc M e // = m — fx^/2a^B 3 . In this picture, the lon- 
gitudinal component of the magnetic field plays the same role as the 
mass does. One has to recall in particular that, in the case of Majo- 
rana neutrinos, the two helicity states do not necessarily have the same 
mass. 



4 A Solvable Model 

We now turn to the description of a solvable model. Two further approxima- 
tions are made in order to accomplish this. First, we consider a model with 
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a reduced number of flavors: we suppress the flavor matrix structure in pre- 
ceding equations altogether. This means that in the case of Dirac neutrinos 
we are dealing with a single flavor of two helicities, one corresponding to an 
active, the other one to a sterile neutrino. In the case of Majorana neutrinos 
we have genuinely two flavors. One hopes that, despite the reduction in the 
number of flavors, the qualitative insight gained by solving the model will 
facilitate the study of more realistic models. Second, we consider the high 
energy limit of the theory developed so far: this is particularly easy in the 
front form of the dynamics. 

We observe that the terms in items 2) and 3) in the preceding Section 
are proportional to e (z — z'). Hence, in a Fourier representation, viz. upon 
writing 

<p (t, z, x A ) = J dkip (t, k, x A ) exp (-ikz) (21) 

and 

= exp (-ikz), (22) 

one finds that at high energies (k ^> m) the terms in items 2) and 3) are 
negligible with respect the local term in item 1). (In the last equation V 
stands for the principal value.) 

Thus, if one is interested in high energy neutrinos, the Hamiltonian can 
be approximated by a local one, given by the first term in (^Of). 

For this reason, it is convenient to work in a coordinate representation 
for the density matrix. The equation of motion then reads: 

— id t (z,x\p(t) z' ,x') = p\Plo ■ B ^x, ^_ j ( z, x\p(t) z', x') 

- pV2(z,x\p(t) \z',x')a-B [x',^S\ (23) 

In this equation, x stands for the transverse part of the coordinate and 
a • B is the two dimensional scalar product in transverse space. Of course, 
the coordinate x 3 had to be expressed by z and t; hence the i-dependence in 
the magnetic field. 

We choose the initial condition so as to describe a neutrino produced at 
x = and with a fixed value of k: 

(z, x\p (0) \z'x') = 5 2 (x) 5* (£) 6XP ik ^~ Z>) Ps (0) , (24) 

where p s (0) is the initial value of the spin density matrix. 

The variable k being large, the function expik(z — z') is rapidly oscil- 
lating unless z ~ z'\ it is permissible to put z = z' in the coefficient of the 
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exponential in (0). Further, in the approximation used, the dynamics de- 
scribed by eq. (p3|) is independent of k and of x. Therefore, the dependence 
of p(t) on these variables is entirely determined by the initial condition and 
the dynamical equation reduces to an equation involving the spin density 
matrix alone, as in non relativistic spin dynamics. From now on, we omit 
the the subscript s and we have: 



id t p (t) = p\Pl 



a-B 



s/2 



P(t) 



(25) 



(Here and in what follows, x = is understood.) 

This equation can be solved by the standard time ordered series, viz. 



p(t) = p(0) 
+ ipV2 



dt' 



a 



B 



V2 



,P(0) 



+ 
+ 



2! 



dt'dt"T 



a-B 



t r 



V2 



a 



B 



V2 



,P(0) 



We choose the initial condition as: 



p(0) = -(l + Sa 3 ), 



(S 2 <1 



(26) 



(27) 



since neutrinos are produced with a definite helicity. (In the case of Dirac 
neutrinos, S = ±1, depending on whether a neutrino or anti neutrino is pro- 
duced. In the case of Majorana neutrinos, S may assume any value between 
the limits stated above, depending on the production mechanism.) 

Next, we average the solution, (|26"D over the magnetic field. We choose 
the generating functional of the moments as follows: 



Cl 



JvBexp- - J d 3 xd 3 x'Bi (x) C^ 1 



x, x') Ba (x' 



exp J d 3 xji (x) Bi (x) 

( x did a 

4tt(B 2 ) 



(L- 2 - V 2 ) 2 6 3 (x-x 1 ) 



[28) 



In the last equation, L and (B 2 ) stand for the correlation length and mean 
square magnetic field, respectively. The measure is normalized such that 
Z [0] = 1. The transverse projector is needed in order to make the correlation 
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functions solenoidal. With the choice of the tensor given in fl28|), the 
leading term in the long distance behavior of the correlation function is oc 
exp — \x — x'\. In order to average equation ( |2T)| ) over the magnetic field, one 
integrates over B 3 and sets the third component of the source equal to zero. 
The transverse generating functional reads: 



T>B exp 



Itt(B 2 ) 



d x 



B\x) ( 8 AB - ] B B 



x exp i / d 3 xj (x) ■ B (x) 



(29) 



We now notice that in the equation (|26|) , terms containing odd powers of 
\i are also odd in B A . Therefore, in the limit j — ► the average of those terms 
vanishes. The even terms in the series are obtained by taking the appropriate 
functional derivatives of (^). All of them are expressed in terms of multiple 
time integrals of Cy (|t — 1'\) and its powers: those integrations are easily 
performed. It is sufficient to illustrate the procedure for the second order 
term in (p6|). 

Carrying out the integrations, one gets: 



-^-S{f dt'dt" 
2 Jo 



a 



B, 



a-B,a 3 ]]) = -fi 2 a 3 (B 2 )tL (l - exp-- 



For large times the result in the last equation is just proportional to t. The 
higher order terms follow a similar pattern. The end result is: 



(P (*)> ~ g i 1 + 5(J 3exp-^ 



(30) 



with 



T 



2fj 2 (B 2 )L. 



(31) 



5 Discussion 



The solvable model described in the last section leads to the remarkable re- 
sult that in the random field the behavior of the helicities is an ergodic one: 
irrespective of what the initial density matrix was, for t > T, the helicities 
are equally distributed. In the case of Dirac neutrinos, this is rather un- 
interesting: roughly 1/2 of them is a sterile one. Unfortunately, however, 
calculated neutrino fluxes emerging from such astrophysical objects as an 
AGN usually cannot be trusted to an accuracy which would permit an ob- 
servational testing of the result. In the case of Majorana neutrinos, however, 
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the two helicity states correspond to two different flavors. Given the fact 
that the neutrinos produced arise mostly from pion decay, the presence of 
the other flavor in roughly equal proportion is an observationally testable 
result. We conjecture that the situation is similar if all neutrino flavors are 
properly taken into account. Should this conjecture be verified by future 
calculations, one would have a very important observational prediction. If 
neutrinos are Majorana particles and they possess substantial (say, 10 _12 /is) 
transition moments, the flavor distribution of neutrino events observed in a 
high energy neutrino telescope would be practically uniform. For this reason 
it would be of utmost importance to 

• observe high energy neutrinos emerging from point sources, 

• develop flavor sensitive detection techniques in neutrino telescopes. 

We remark that the question of flavor change in the case of magnetized 
Majorana neutrinos can be investigated in terrestrial neutrino experiments 
as well, see ||. The two types of experiments referred to above are comple- 
mentary to each other. If both types indicate the presence of flavor flipping, 
one would have a window on post standard model physics as well as the 
structure of important point sources of high energy neutrinos. 

From the theoretical point of view, the main result of this work is the 
development of a consistent formalism in order to treat the propagation of 
high energy neutrinos in a random magnetic field (more generally, in a ran- 
dom environment). This problem has been treated before, see for instance, 
ref. and the virtually complete bibliography quoted there. (The authors 
quoted in those references concentrate mostly on the early universe and the 
role of sterile neutrinos in it[] The main advantage of the formalism devel- 
oped in refs. |3|, [| and adapted to highly relativistic systems in this paper is 
twofold: it is guaranteed to be free of internal inconsistencies and it allows 
the development of controllable approximation methods. 

A final comment is in order. Assuming that one observes flavor conversion 
in some AGN, one may ask whether the theory developed here can be used 
to test some properties of the magnetic fields in the source. We argue that 
the answer is in the affirmative to some extent. The measurable quantity 
characterizing the magnetic field is the coefficient given by equation (pTj) . 
One notices that this expression contains the quantities characterizing the 

4 In most of the works quoted, an average is taken at the level of the evolution equations 
rather than their solutions. It has been known for some time that averaging the evolution 
equations over a random variable is an incomplete procedure (important physics may be 
missed). It may also lead to internal contradictions. For a lucid exposition, see JX0| - A 
notable exception is ref. jy], where a simplified equation is given a correct treatment. 
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magnetic field only in the combination fi 2 (B) 2 L. This is not an accident. 
There is an important class of distributions of the magnetic field characterized 
by the following: 

1. The magnetic field has a zero mean value, 

2. The distribution is approximately Gaussian (higher cumulants are ap- 
proximately zero), 

3. the distribution contains only one length scale. 

Every such distribution leads to a damping coefficient of a form which can 
differ from ( |31~]) only in a numerical factor. In essence, this follows from 
dimensional analysis. 

In fact, the magnetic moment and the magnetic field strength enter the 
Hamiltonian in the combination fiB; since (B) = 0, this combination must 
enter a physical observable quadratically. Finally, in order to form a quantity 
of inverse length (the damping coefficient), the expression has to be multiplied 
by a characteristic length in the distribution. 

As another example, Enquist et. al. , [0 derived a damping coefficient 



of this form, even though they use a field distribution quite different form 
oursQ. (The distribution in |T^] has a vanishing correlation length, but a finite 



domain size.) Thus measuring such quantities as the damping coefficient can 
only yield information about the class of the magnetic field distribution, not 
about its details. 
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